gusucode.com > 《MATLAB图像与视频处理实用案例详解》代码 > 《MATLAB图像与视频处理实用案例详解》代码/第 19 章 基于语音识别的信号灯图像模拟控制技术/melbankm.m
function [x,mn,mx]=melbankm(p,n,fs,fl,fh,w) if nargin < 6 w='tz'; if nargin < 5 fh=0.5; if nargin < 4 fl=0; end end end f0=700/fs; fn2=floor(n/2); lr=log((f0+fh)/(f0+fl))/(p+1); bl=n*((f0+fl)*exp([0 1 p p+1]*lr)-f0); b2=ceil(bl(2)); b3=floor(bl(3)); if any(w=='y') pf=log((f0+(b2:b3)/n)/(f0+fl))/lr; fp=floor(pf); r=[ones(1,b2) fp fp+1 p*ones(1,fn2-b3)]; c=[1:b3+1 b2+1:fn2+1]; v=2*[0.5 ones(1,b2-1) 1-pf+fp pf-fp ones(1,fn2-b3-1) 0.5]; mn=1; mx=fn2+1; else b1=floor(bl(1))+1; b4=min(fn2,ceil(bl(4)))-1; pf=log((f0+(b1:b4)/n)/(f0+fl))/lr; fp=floor(pf); pm=pf-fp; k2=b2-b1+1; k3=b3-b1+1; k4=b4-b1+1; r=[fp(k2:k4) 1+fp(1:k3)]; c=[k2:k4 1:k3]; v=2*[1-pm(k2:k4) pm(1:k3)]; mn=b1+1; mx=b4+1; end if any(w=='n') v=1-cos(v*pi/2); elseif any(w=='m') v=1-0.92/1.08*cos(v*pi/2); end if nargout > 1 x=sparse(r,c,v); else x=sparse(r,c+mn-1,v,p,1+fn2); end